Effective dynamics of a nonabelian plasma out of equilibrium 
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Starting from kinetic theory, we obtain a nonlinear dissipative formalism describing the nonequi- 
librium evolution of scalar colored particles coupled selfconsistently to nonabelian classical gauge 
fields. The link between the one-particle distribution function of the kinetic description and the 
variables of the effective theory is determined by extremizing the entropy production. This method 
does not rely on the usual gradient expansion in fluid dynamic variables, and therefore the result- 
ing effective theory can handle situations where these gradients (and hence the momentum-space 
anisotropies) are expected to be large. The formalism presented here, being computationally less 
demanding than kinetic theory, may be useful as a simplified model of the dynamics of color fields 
during the early stages of heavy ion collisions and in phenomena related to parton energy loss. 
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I. INTRODUCTION 

The results of numerous experiments on ultrarel- 
ativistic heavy ion collisions carried out at the Rel- 
ativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) clearly point to the conclu- 
sion that a hot, dense and opaque nuclear medium is 
created in such events [lHSj- By now, the standard 
picture for the evolution of such matter consists in 
a pre-equilibrium stage dominated by strong chro- 
moelectromagnetic fields, followed by the formation 
of the Quark Gluon Plasma (QGP) in local thermal 
equilibrium. The QGP thus formed expands and 
cools under its own pressure, going through the de- 
confinement transition in which hadrons are formed 
and later on detected (for reviews see 0-Q and ref- 
erences therein). 

Before the collision, the occupation number of glu- 
ons at high energy is so large that an approximation 
in terms of classical gauge fields obeying Yang-Mills 
equations becomes more suitable than a description 
in terms of on-shell particles 0, [l(| ■ The time scale 
for this highly nonlinear regime occurring at the ear- 
liest stage of a heavy ion collision is Qj 1 ~ 0.2 fm/c 
at RHIC, where Q s is the saturation scale. 

The dense system appearing between first impact 
and the formation of the equilibrated QGP develops 
Chromo-Weibel instabilities (see e.g. [ll|, [Hi): the 
non-equilibrium anisotropic distribution of the par- 
tons is responsible for the fast growth of the chro- 
momagnetic plasma modes, which in turn isotropize 
the system and speed up the thermalization process 
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to yield a thermalization time ~ 1 fm/c. This is 
roughly the value of the thermalization time that 
can be inferred from a comparison of state-of-the- 
art hydrodynamic simulations to data. Moreover, a 
fast parton transversing the matter formed in heavy 
ion collisions excites color fields and loses energy, a 
phenomenon known as jet quenching (l3l Il4j |. For 
these reasons, the dynamics of classical gauge fields 
(both coupled or not to hard partons) has received 
a lot of attention in the context of heavy ion phe- 
nomenology 

The color fields "see" the evolving matter (which 
in our case is composed of colored partons) through 
the conserved color current, as dictated by Yang- 
Mills equations. The dynamics of the conserved 
currents, i.e. the color current and the total energy- 
momentum tensor of the combined system of matter 
plus gauge fields, is in general very complex, since in 
principle it must be computed from the dynamics of 
the matter fields (or particles in the kinetic limit). 

In the semiclassical kinetic approach, the coupling 
of hard and soft modes is implemented through a 
nonabelian transport equation [5o1 . [56j which deter- 
mines the evolution of a one-particle distribution 
function f, which is a matrix in color space (see 
Sectio n III B I) , or equivalently through Wong's equa- 
tions [57|; see also Refs. |58r-69]. Once f is known, 
the color current acting as the source in Yang-Mills 
equation is completely determined, so the dynamics 
of the gauge fields can be found. 

Here, we take the view that it is a reasonable hy- 
pothesis to assume that the dynamics of the con- 
served currents is largely determined by the conser- 
vation laws themselves. This opens up the possibil- 
ity of constructing an effective theory incorporating 
the conservation laws, that allows one to investigate 
relevant aspects of the dynamics of gauge fields (for 
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example, the border between stability and instabil- 
ity, or the back reaction of hard particles on the 
evolution of color fields) in a simpler context as com- 
pared to the microscopic theory. 

The degrees of freedom of the effective theory in- 
clude hydrodynamic variables such as flow velocity. 
However, since we are working in a strongly out of 
equilibrium regime, it is not possible to obtain a 
closed dynamics based on hydrodynamic variables 
alone. 

We obtain a closed theory by identifying a ten- 
sor X^ (which is introduced in Section UTTl) through 
which the system couples to the hydrodynamic vari- 
ables. One can think of the introduction of the 
nonhydrodynamic variable A Ml/ as a simple way of 
modeling the back reaction of the distribution func- 
tion (which could correspond to a highly nonequi- 
librium situation) on the hydrodynamic modes (i.e. 
those modes associated to conservation laws and 
hence relaxing much more slowly) . A somewhat sim- 
ilar situation is discussed in |70M74j in the context 
of the Entropy Maximum Principle (EMP), in [75| 
(see also [76|) where the moments of the collision 
term of the classical Boltzmann equation are inter- 
preted as independent variables rather than as infi- 
nite moment series, in [77] and [78| in the context of 
"anisotropic hydrodynamics" , in |79f in the frame- 
work of divergence-type theories [8Q| (see also [8ll — 
HIj])j m 01 where the closure is obtained from an 
expansion of the distribution function in moments at 
all orders, and in [85l |86| in relation to the Entropy 
Production Principle (EPP). 

Here, the tensor X^ is identified with a La- 
grange multiplier in a well-defined variational prob- 
lem whose solution yields the distribution function 
that extremizes the entropy production (for a review 
of this method see [87| )• This procedure results in a 
closure for the distribution function of colored par- 
ticles which has two satisfying properties: (i) it is 
nonlinear in the variables of the effective theor y, th us 
generalizing Grad's quadratic ansatz [H, l89l |94{ in 
a nontrivial way [H, [9(| , and (ii) it does not rely in 
any way on the gradient expansion in fluid variables. 
The equation of motion for X^ v is then obtained from 
this closure by the method of moments. 

The result is an effective theory for the dynamics 
of color fields coupled to colored particles that can 
handle highly nonequilibrium situations, for exam- 
ple the large momentum anisotropy present at early 
times in heavy ion collisions. Our formalism is a 
simplified model that shares with the true dynamics 
the conservation laws, Lorentz and gauge invariance, 
being causal, and satisfying the Second Law. The 



limitations of our model will be discussed later on. 

We would like to stress that our approach is not 
new. The idea of obtaining a closure from a varia- 
tional principle whose Lagrange multipliers are iden- 
tified with macroscopic variables is at the heart 
of the so-called Extended Thermodynamic Theories 
(see e.g. [Hj]), which use as a variational principle 
the EMP. However, in most applications of this the- 
ory the EMP is only used to express / in terms of the 
Lagrange multipliers but not to obtain their dynam- 
ics. In Refs. fmW^ the EMP was used to accom- 
plish both tasks in the context of electron transport 
through mesoscopic semiconductors in the nonlinear 
and nonequilibrium regime (i.e. under conditions of 
very strong electric fields and large gradients). The 
result was a closed effective theory which could re- 
produce well the results of kinetic Monte Carlo sim- 
ulations. 

In this paper, instead of using the EMP we rely 
on the EPP, which allows us to find the distribu- 
tion function that extremizes the entropy production 
subject to constraints on the conserved currents pS7j . 
We note that the EPP was used in [85[ to construct 
a model of nonequilibrium electron transport and in 
(86j to describe radiative heat transport in a pho- 
ton gas, obtaining results that agree well with more 
sophisticated approaches. 

The motivation for the choice of the EPP is three- 
fold. This principle is deeply connected to nonequi- 
librium statistical physics, in particular to the evo- 
lution of fluctuations around a stationary state (we 
shall not discuss this connection here; for details 
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931V Moreover, as it is rigorously 
87] for nonrelativistic systems (see 
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shown in Ref. 

also [Isl), the distribution function obtained from 
the EPP actually solves the linearized Boltzmann 
equation, which means that the closure provided by 
this variational principle is able, at least in princi- 
ple, to capture some features of the microscopic dy- 
namics. Finally, the transport coefficients obtained 
from the EMP can differ substantially from those 
computed from the EPP, the latter providing bet- 
ter agreement with the results obtained from kinetic 
theory (see for instance [zO-tZl H, H, H3 ) ■ 

In the spirit of finding the simplest possible theory 
incorporating the conservation laws, we postulate a 
simple BGK form for the collision term (94| in the 
Boltzmann equation, and obtain solutions at second 
order in the relaxation time. By going to second 
order in the relaxation time (instead of the usual first 
order treatment) we expect to broaden even more 
the range of applicability of the effective theory, to 
be able to describe better those situations with large 
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momentum-space anisotropics. 

In relation to our developments, we note in par- 
ticular Ref. [1(| in which the so-called "anisotropic 
hydrodynamics" (that was developed in [zJlzl]) is 
coupled to color fields in a way reminiscent of mag- 
netohydrodynamics. In j36| , an effective model for 
the dynamics of color fields coupled to particles is 
obtained from the Boltzmann-Vlasov equations by 
the method of moments and then applied to study 
chromoelectric oscillations in a dynamically evolving 
anisotropic background; see also [35j . As it will be- 
come clear in what follows, our approach resembles 
those adopted in these studies in that they are effec- 
tive theories capable of dealing with large deviations 
from equilibrium. 

The colorless version of the effective theory devel- 
oped in this paper was investigated in [8l[ (although 
there it was obtained in a different way as the one 
followed here). It was then applied to study the 
evolution of matter created in heavy ion collisions 
at RHIC and to the calculation of hadronic observ- 
ables in [Hf, and compared to second order fluid 
dynamics, to which the colorless effective theory re- 
duces when deviations from equilibrium are small. 
The connection between the (colorless) effective for- 
malism and kinetic theory was established in 83] . 

Similarly to what happens in the colorless ver- 
sion of the effective theory [gj [Hf, the formal- 
ism presented here reduces to the so-called "chro- 
mohydrodynamics" [13, [H, H(| HU when deviations 
from equilibrium are small (see Section |IV|) . Previ- 
ous studies based on chromohydrodynamics include 
the calculation of the wake potential induced by a 
fast parton [1^ as well as collective excitations and 
instabilities in the QGP dJii IM IhI ■ 

This paper is organized as follows. In Section 
lrfl we describe the basic theoretical setup for our 
developments, give a very brief overview of the ki- 
netic theory of a nonabelian plasma, and introduce 
the conserved and entropy currents. In Section Mil 
we obtain a closure for the distribution function 
from the entropy production principle and derive 
the evolution equations of the effective theory by the 
method of moments applied to the transport equa- 
tion; this Section contains our main results. In Sec- 
tion |IV] we compare our developments to the matrix 
approach to chromohydrodynamics, and as a sim- 
ple illustrative example we compute the polariza- 
tion tensor of the colored plasma including a finite 
relaxation time for the fluctuations. We conclude 
in Section [V] with some comments on the possible 
application of the developed formalism to the dy- 
namics of color fields and instabilities in heavy ion 



collisions. 



II. THEORETICAL SETUP 
A. The system 

We are interested in obtaining an effective theory 
for the dynamics of a system of colored particles in- 
teracting with nonabelian classical gauge fields. In 
this work we will deal with scalar particles coming in 
three colors, which in our simple model would rep- 
resent massless and spinless quarks. We shall there- 
fore consider a classical Yang-Mills field coupled to 
conformal scalar matter in the fundamental repre- 
sentation of SU(3). 

In what follows, we will use (/i, v,...) to denote 
world indices and (a, 6, . . .) to denote internal (color) 
indices. We shall denote with N — 3 the dimen- 
sion of the fundamental representation, and use n 
to indicate a generic dimension (n = 3 or n = 8 for 
the fundamental or adjoint representations, respec- 
tively). The generators T a are traceless hermitian 
n x n matrices with commutation relations 

[T a ,T b }=iC c ab T c (1) 

and trace 

trT a T fc = hSab (2) 
The Yang-Mills field is A M = A£T a . The field tensor 

F M „ = d^A,, - d y A^ - ig [A M , A v ] (3) 

belongs to the adjoint representation of the gauge 
group. The equations of motion for the Yang-Mills 
field are 

D M F^ = -Q [J,] (4) 
where the covariant derivative is 

D M X = S p X-i 5 [A M ,X] (5) 
Q is a projection operator 

Q [X] = 2 V T o trT X = X -trX (6) 

L — ' n 
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Eq. ((D) implies the Bianchi identity 

Q [D M J"] = (7) 
We also have the energy-momentum tensor 

rpHU _ . rriflV / o\ 

YM — Y M \°) 

where 



T^ M = F^F^ - -g^F x PFx P (9) 

is traceless in world indices. Using the identity 
D( M F„,\) = (where brackets mean symmetriza- 
tion) we get 



C. Conserved currents 

We now introduce the matter entropy and the con- 
served currents of the microscopic theory. The non- 
abelian current reads 

Ja =9 J Dpp x f (13) 

where Dp = d 4 pS (p 2 ) / (2tt) 3 . 

We have not written down an explicit equation 
for the matter stress-energy tensor T^", but we know 
that since the total stress-energy must be conserved, 
we must have 

T% v = -T^ M . u =trJ X F>* X (14) 

We get this by writing 

ir = trir (is) 



tr{D M T^ M +J A F- A } =0 



(10) 



B. Kinetic theory 

In principle the scalar matter should be described 
by quantum field theory. The reduction of the 
nonequilibrium quantum field description to kinetic 
theory is fairly established by now, see e.g. [H| and 
references therein. 

The kinetic equation which governs the evolution 
of the one-particle distribution matrix f reads [EH 
HI (see also [HHP 



D M f 



df Of 

Fui/"^ "1" ~7{ F uis 

ap v op v 



where 



D M f 



= sign I co i 
(11) 

(12) 



with expressed in the fundamental representa- 
tion, f (X,p) is an N x N matrix (N — 3 for quarks) 
and obeys ft (X,p) = f{X,p). 

The collision kernel on the right hand side encodes 
the interaction among the hard excitations of the 
matter field, and is, in general, a complicated func- 
tional of the self-energy [58]. However, to carry out 
our developments we really do not need to consider 
it in much detail, and we anticipate that later on we 
will use a phenomenological linear collision operator 
that will suffice for our present purposes. 



Dp p^p w { 



(16) 



In Eq. ([14]) , the semicolon stands for an ordinary 
derivative. We shall drop the subindex m for 
and write T^" in what follows. 

Eqs. ([7]) and (jTUJ) are identically satisfied provided 



J Dp sign (p°) tr(T a I coi ) = 
tr J Dp sign (p a )p^I col =Q 
The entropy current is 



(17) 



S" = J £>p^sign(p°)tr{(l + f) ln(l + f ) - f lnf} 

(18) 

leading to the entropy production 

/ Uptr^lnf-^l + f)} (19) 



We note that to go from Eq. (|T8} to Eq. dT9J> 
one must assume that [f, D M f] = (see HH). If this 
condition is not imposed on the distribution func- 
tion, the entropy production contains terms of the 
form tr{F^ t „, f } which contribute to entropy produc- 
tion even in mean field. Assuming [f , D M f] = then 
corresponds to assuming that entropy is produced 
solely by collisions among the particles. We shall 
stick to this approximation in what follows. 
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III. EFFECTIVE THEORY 



where we have used that 



A. Obtaining a closure for f 

We will now obtain an expression for the one- 
particle distribution matrix in terms of variables of 
the effective theory. To do so, we will rely on the 
EPP [87] discussed in the Introduction. 

This method allows one to find the distribution 
function which extremizes the entropy production 
Sft, subject to the constraints that the conserved 
currents take on known values. It provides a pre- 
scription to associate a distribution function to given 
macroscopic currents, yielding a nonlinear closure 
that generalizes the well-known Grad's quadratic 
ansatz [H, [H, HI [23 in a nontrivial way [H, Hoj]. 
The EPP does not rely on a gradient expansion, 
which is usually invoked when deriving hydrodynam- 
ics from kinetic theory. This results in effective theo- 
ries capable of describing highly nonequilibrium and 
nonlinear situations quite reliably as compared to 
microscopic approaches [7(il-l74l HH-HtJ . 



J Dp tr {I co , m(/o)} = 



(25) 



Given the expression for the entropy production 
given by Eq. (|2~4l and taking into account that we 
will linearize the collision operator, in order to sat- 
isfy the H-theorem we introduce a new variable Z 
such that 

e z = l + X (26) 
We will work to quadratic order in Z, so we get 



Parametrizing Z as 



z = — c (0) i + C a T a 



(27) 



(28) 



and using that 



1. Deviations from the unperturbed state 

For simplicity, in what follows we will neglect 
quantum statistics. Given tr(T A " y ) = T Miy , we can 
define a flow velocity w M and a temperature T by 
using the Landau-Lifshitz prescription. We have 



(20) 



where p{T) is the energy density as obtained from 
the equation of state. The prescription amounts to 
matching the local nonequilibrium state of the flow- 
ing real matter to a fiducial perfect fluid. We will 
show later that the stress tensor II MI/ is transverse, 
i.e. w At IP" y — 0, so the prescription is consistent. 
With and T, we can define /3 M = u^/T and then 
construct 



fo(x»,P 

We write the distribution function as 



(21) 



f = / [l + (l + /o)xWo[l + x] (22) 
The entropy production 



(23) 



then reads 



S u* = -f ^Ptr{I co; lnf} 
S*=-f £>ptr{I coi ln(l + X )} (24) 



with 



K c b = ^(iCabc + dabc) 



(29) 



(30) 



where d a b c are the symmetric structure constants, 
we get 

z " = ^(C (0) ) 2 i+^C (0) C a T Q +^C Q C Q i+^C a C b df t c 

(31) 

The quantities £ a can be identified with color fugac- 
ities ( a = fJ. a /T, where /i a are the color chemical 
potentials needed to conserve color (see Eq. ([T7)) ). 
Having Z, we obtain 



y= jL((°)i + ^T a + i 
s 2 



— (C (0) ) 2 i 



2 1 1 

___ /"(0) /-0,rr\a _|_ C a C a \ \ — -C a d CL ^ r J? C 

N 2N 2 c 



(32) 



2. Currents, entropy production and collision term 

To solve the variational problem we must express 
the currents and the entropy production in terms of 



G 



The shear tensor 
reads 



(33) 



(34) 



where we have introduced the notation 

Dpf (...) = <(...)) (35) 



From Eq. ([13]) we get the expression for the color 
currents 



1 



where we have defined n = g (u>) . 
For the entropy production, we have 



(36) 



St 



Dp tr(I coi Z) 



(37) 



The collision operator contains color independent 
and dependent parts, so, similarly to the decompo- 
sition used for Z, we put 



T 7 — r^l -I- T a T Q 



(38) 



We then have 



tr(I coi Z) = ^C (0) + i^C a (39) 
where we have used Eq. ([2]), so 

^ = - J Dp ^C (0) + ^ oi C Q ) (40) 



We shall write the collision operator as 
Ip = ~R[FR[Z[] 



where we have put 



Icoi = /ol^ 



(42) 



In Eq. (l4"Tj) , r is the relaxation time, F — F(u) is 
an arbitrary function of energy, and R is an operator 
enforcing the integrability conditions given in Eqs. 
([TT]) . The projector R explicitly reads 

R[Z] = Z-R ab (tr(T„Z)T 6 ) - 1R„ (p"tr(Z)) (43) 



with 



1 



1 



JVM"" ' iVM 



and i?ab = 26 a b. Here, 



(44) 



(45) 



is the spatial projector. Note that this form for 
Ip guarantees that the Second Law holds exactly. 
Moreover, it is flexible enough to include the im- 
portant cases of Marie's relativistic generalization 
of the BGK model, corresponding to F(w) — T, as 
well as the Anderson- Witting model corresponding 
to F(ui) = lo (for a discussion of these models in con- 
nection to the freeze out stage in heavy ion collisions 
see jsil23|). 

From Eq. (14"TT) we then have 



and 



with 



and 



r(o) 



T" 



1 



2riV 



R[FR[C 



(0)1 



-—R[FR[C}] 



R[C (0) ] = 1_ C (0)_ R ^ ^ C (0)^ 



1 



TaCa ~ ^ Tc (Cb) 



TaR[(a 

The integrability conditions then read 



= 



and 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



(41) Note that there are no restrictions on (p p Ig\ or on 



(0) 



3. Variational equations 

In our case the conserved currents are J afJ/ and 
T^ v , so that the variational problem becomes 



SZ 



" V T/1 * " V J* M ] = (52) 
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where and A aM are Lagrange multipliers forcing 
the energy momentum tensor and the color currents 
to take on their known values. 
We then have 



h A., ,, 



— Ajui/— ptj; h A a 



' <5C (0) 



8J a » 



(53) 



Using Eqs. flM|), © and together with Eqs. 
(|4*o) and (|4"T|) . the variational equations (|53")) become 



1 



R[FR[C W }] = A 



(o) 



and 



1 

N 
(54) 



-R{FR{( b }} = A^pV^+sW 



' (55) 

The solution to the variational problem entails 
finding Z as a function of the Lagrange multipli- 
ers, and the outcome is a closure for the distribution 
function f, i.e., an expression for f in terms of the 
variables of the effective theory. 

We will solve the equations (|5"4"]) and ([531) to second 
order in the relaxation time r. To this end, we will 
follow [83| and expand the nonequilibrium correction 
Z and the Lagrange multipliers as follows 

C (0) = C (0) + C (0) 

j~a /-a i /-a 

S — SI T S2 

A -A (0) +A (1) 

K = K {0) + K {1) 



(56) 



We now go over to solve the equations at first and 
second order in r. We note that the solution given 
here closely follows the one given in [83]. For the 
reader's convenience, a brief summary of the logical 
steps carried out to obtain the closure for the distri- 
bution function can be found in Section (IIII A6|) . 



4- First order solution 
At first order, the variational equations are 

-Lr[FR[^}] = AWp'V (57) 

and 



1 



R[FR[£]] = 



(58) 



Eqs. (|4*^|) and (|58"j) imply that £i must be indepen- 
dent of p u , as it must be given that £i are color 
fugacities. 

From Eq. (1571) we get 



<pW> = 



(59) 



Setting p — k we obtain — 0. Without loss of 

generality we can take A$ = 0. Since we are dealing 
with a conformal theory (p l p>) = S 1 ^ (w 2 ) /3, and we 

get a|°^ = 0. We thus obtain 



(60) 



The first order shear tensor is then 

nr = (pVd 0) ) (6i) 

We find that n 00 = II fe0 = 0. Using that 

(G(w)pV'pV) = 1 (6 %1 S kl +S %k S jl +S u S 3k ) 

15 

for any function of energy G(u>), we get 
2tN /lj 4 ' 



Hi 



y _ 



15 \ F 



_ \ A (o)« 



(63) 



which is traceless and transverse. The first order 
color current is (recall that must be independent 
oip u ) 



■7?" = </ 



(64) 



The lowest order nontrivial contributions to the 
entropy current and the entropy production are 

(65) 

in the rest frame, and 



respectively. The entropy flux (1651) and the entropy 
production (156")) computed from the first order so- 
lutions (Ci°\Ci) to the variational equations are al- 
ready quadratic in deviations from equilibrium, so 
we do not need to consider higher order contribu- 
tions. Note that, because of the first order equation 
(|58|) . £i does not contribute to the entropy produc- 
tion at quadratic order. 
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5. Second order solution 
The variational equations at second order read 

-Lr[FR[^]} = Ag,W + ^y^cf 



and 



AT 



N \ P (,1 



(67) 



±R[FR[$]] = g^\ b ^p^i 0) + X^P^Ct (68) 
The integrability conditions then read 

rAA« (pW+l ((d 0) ) 2 Fo;)+ g rA o(0) C? <- 2 ) = 

(69) 

riVA« (pV'o;) + ffrAf <pV) = (70) 
corresponding to Eq. (|50p and 

^Ar(^ O) )=0 (71) 

corresponding to Eq. (|5T]l. The latter equation 
shows that A^ " 1 = 0. 

Any term in A^ 1 ,] which is not strictly required 
by the integrability conditions can be absorbed into 
Ajui , so there is no loss of generality if we can take 



A 



(i) 
oo 



0. Moreover, from Eq. (|70|) we see that A 



and A 



a(0) 



- is a solution to the integrability 
condition. Therefore, we can write 



where we have put 



A 



1 



A 2 t (cj3) 
The quadratic equations then read 



(72) 



(73) 



i?[Fi?[cf ]] = -tJV Acj 2 + ^(CD^ (74) 



1 



integrated agaist cj/F. Therefore, the solution to 
Eq. (d) is 



C (o) = l (c (o) )2 _ 7VrA ^!_^ 
S2 N F F 



where 



.4 



NtA 



Similarly, we obtain from Eq. (1751) 

?6 



<2 



where 



I?" 



/ ,^ DO 



(76) 



(77) 



(78) 



(79) 



Note that A oc (d **) 2 is already quadratic in r, 
which means that £5 is actually third order and 
therefore can be neglected. 

We are ready to compute Tlrf and J^; see Eqs. 
dMD and HMD- We have 



X1 2 



2A Ul j 4^^ 



Note that n° ,; 



NrK^ - A— 
F F 

(80) 



0, but U9, ^ 0, so is not 
the true correction to the energy-momentum ten- 
sor, whereby the parameter T in our equations is 
not the true temperature (that would be measured 
by an observed moving with the local rest frame). 
To obtain the physical correction to T^ v , which we 
call ^2phys' we mus t perform a temperature shift. 
Putting T = T phys - ST with 8Tdp/dT = Hf\ the 
true correction reads 



n 



2,phy; 



= 1$ - -8 i m 



1)1] 



(81) 



which is traceless. Taking into account the above we 
find 



11 2,phys 



3 lml (82) 

+ T(u) A^CKi 



and 



with 



R[FR[( b 2 }} = -tAu; 2 ^ 



(75) 



Note that the left hand side of Eq. (|74p vanishes 
when integrated against to (which means that the 
equation has a solution) but does not vanish when 



L 



T 2 N /uj 6 
~A5~\F^ 



(83) 



To obtain the last term in Eq. (|82[) we have used 
the identity (lu 2 ) = 2T(w). Note also that Tl 2phys 
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is transverse. For simplicity, in what follows we will 
drop the subindex "phys" in lC° 



2,phys ' 



The quadratic contribution to the color current 
can be computed directly from Eqs. (1551) . ([75]) and 
([75)1 . We get 



(84) 



where we have used that X^-' — so that the last 
term in Eq. (l36l) drops out. 

Using the closure picked out by the EPP, we 
have completed the task of expressing the currents 
tr(T MI/ ) and J M in terms of the variables of the ef- 
fective theory. 



6. Summary of the EPP method 

For clarity, we now briefly summarize the main 
logical steps followed to obtain the closure for f . 

From a linear transport equation, we set up the 
variational problem given in Eqs. (|53|) . The solu- 
tion to these equations gives the distribution func- 
tion that extremizes the production of entropy sub- 
ject to the constraints that the conserved currents 
take on known values. Using the expressions for the 
stress tensor, the color currents and the entropy pro- 
duction given in Eqs. and (|4"0")l . respec- 
tively, together with the integrability conditions for 
the collision term, Eqs. PS)) and (14"T1) . the varia- 
tional equations become Eqs. (|54"|) and ([55]). 

We then expand the Lagrange multipliers and the 
nonequilibrium correction Z in powers of the relax- 
ation time t (Eq. ([56]) ). and solve the variational 
equations to second order in r. The result is an 
expression for the correction Z in terms of the La- 
grange multipliers, which is given by Eqs. (|60[) and 
(|76p . From these equations, we can express the shear 
tensor and the color currents in terms of the La- 
grange multipliers, obtaining Eqs. (1521) 
and f84]). 

The distribution function then reads 



f = /o 

N 2 t 2 



1 + 



2F 2 
+ /o 



(85) 



have been used in different contexts to obtain clo- 
sures for the distribution function, which resulted in 
effective models whose dynamics compared well with 
kinetic theory (see e.g. I87l . 192} for a broad perspec- 
tive and 

[Z^-EE HI, EE [86] for specific applications). 
To obtain a dynamical model, we must now find 
the equation of motion of the Lagrange multipliers. 



B. Equations of motion 

The EPP described above has provided us with an 
expression for the distribution function f in terms of 
the reduced set of variables of the effective theory: 
the usual hydrodynamic variables u M , p, the nonhy- 
drodynamic tensor Ajfy encoding the back reaction 
of f on the hydrodynamic modes, and the color fu- 
gacities Ci ■ This expression is given in Eq. (|55|) . 

The equations of motion for u M , p and are 
the conservation equations for T^ v and J aA1 , re- 
spectively. However, the conservation laws are not 
enough to fully determine the dynamics of the sys- 
tem, and an equation governing the evolution of 
Ii^ u must be given. Usually, this is done within 
the grad ient expansion for hydrodynamic variables 
[7l.l8ll88l.l89j. Here, instead, li^ is an algebraic func- 
tion of X^ ij and d a , so we must obtain the evolution 
equation for . We will get the latter from the 
kinetic equation. 

It will prove convenient to express our results in 
terms of a new relaxation time (related by a con- 
stant to the previously introduced r) 

(86) 



f 



(i 



ui 
F 2 



and the shear viscosity 

_ 6r n N I uj 5 
V ~ 15T XF 2 



F 



(87) 



which naturally arise in the context of second order 
fluid dynamics 7, 8, 98, 99]. For simplicity, we also 
introduce a new variable 

= tjTX^ 11 (88) 

We shall deal with the conservation equations 
first, and then go over to discuss the evolution equa- 
tion for 7 l - J . 



Conservation equations 



As discussed in the Introduction, this method, or 
very similar ones based on maximizing the entropy, 



Inserting the expression for f given by Eq. ([85 
into Eqs. (0 and (10]), and using that Tg? v = 



10 



—Ty J M v — trJAF MA in the latter, we get equations 
of motion for the velocity u M , the energy density p 
and Ci ■ 

For the matter energy-momentum tensor T^ u = 



1 



n^, where 



J- n 



1 



is the perfect-fluid energy-momentum tensor, 



nf = 



(89) 



(90) 



and 



nf 



rj 2 T 2 
- T (w) A^Ci Ci 



3 7 7(Tp y (9i) 



we get 



j O;i/ + u 1;i/" 



n 



2;i/ 



:«(Ci e + 7 ^Ci c Ci d )«AF^ A (92) 



To avoid being cumbersome, we shall not write down 
the explicit expression for T.^f, but it follows imme- 
diately from Eqs. (l89l)- (pTj) .' 

The conservation equations for u M and p as given 
in Eq. (192|) involve the time derivative of and 
(these quantities appear in the LHS of Eq. (p?2"j) ). 
The evolution equation for is obtained from the 
(covariant) conservation of J M . We get 



(93) 



2. Evolution equation 

We shall now deal with the evolution equation for 
the nonequilibrium tensor 7^. We will obtain this 
equation from the second moment of the singlet sec- 
tor of the transport equation Eq. (|11[) . i.e. from the 
equation that results from taking the trace of (fTTj) . 
The reason for considering the colorless part of the 
transport equation will be discussed after presenting 
the evolution equation. 

The second moment of the colorless part of the 
kinetic equation reads 



Dp p p p a p tl 



tr(D^f)- 2tr F 



df df , 

dp„ dpj, 



Dp sign(p )p> CT tr(I c 



(94) 



We have that p M tr(D Jtl f) = p M 9 At tr(f), where 

tr(f) = f = f (N+cl 0) +Cf+^(Ci (0) ) 2 +|Ci a Ci a ) (95) 

is the colorless part of f at quadratic order. Explic- 
itly, we have (see Eq. ([85])) 



f = fo 
N 2 t 2 



1 + 



, p2 pWp^^ + ^Q 



(96) 



The second term in the left-hand side of Eq. (IM|) 
becomes 



(97) 



where 



/ a = /o(C? + jAgW <? + i ^/Ci e d / ) (98) 

is the color part of f at quadratic order (see Eq. 
\). The right-hand side of Eq. ([M]) is 



N J Dpp p p <T sign(p°)I, (0) 



(99) 



Using the explicit expressions for / and f a given 
in Eqs. (IM|) and in the moment equation (IM|) 
we obtain (recall that ^ = TjTX^- ' 1 ^ and that, for 
simplicity, we assume that 77 is a constant) 



Tim J Tij)/* 



T 

iV 
2^T 



jpa /-a 

2rjT 



ijlinpa 



(100) 



For brevity, we have defined the following quantities 



with 



— — — */■ -Y; fjf> a ^ lm 



(101) 
(102) 
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In order to find an explicit evolution equation for 
A;??„ , we must be able to invert the tensor with which 
it is contracted, namely 

TjT 

For our present purposes it is enough to display the 
equation of motion for the nonequilibrium tensor j lJ 
to linear order. 

Using suitable generalizations of Eq. (frJ2|) to com- 
pute the N's explicitly we get (in the local rest 
frame) 

r = -aid a A'% a - .1(1 + -xKty j 

T " (104) 

+ a 2 Ci a ^ (l 7 Qj) 1 13 

In Eq. (|104[) , a p<T is the first order shear tensor 

a pa = V <p u a> (105) 

where C <pv> denotes taking the traceless and trans- 
verse part of a tensor C and V M = A Miy i9„ is the 
spatial gradient. We have denoted the convective 
derivative by an overdot, i.e. C = u a d a C. The 
parenthesis around indices denote symmetrization. 

The transport coefficients that appear in the 
evolution equation are 

The transport coefficients a\ and ai are novel coeffi- 
cients that couple the nonequilibrium tensor 7 pcr to 
color degrees of freedom. A nonvanishing a\ implies 
that the varying color chemical potentials affect the 
evolution of 7 P<T . The term containing a 2 represents 
the coupling of 7 pcr to the gauge fields. 

If 7V — >■ and Ci = in Eq. (11041) , we recover the 
(colorless) Navier-Stokes limit with r y plJ —t —rja fJ,L '. 
On the other hand, we have already shown in [8l| 
that if we expand 7 P " to second order in velocity 
gradients, the above formalism (with = 0) goes 
over to the second order conformal hydrodynamics 
that was derived in Refs. (HHl]. 

The equations of motion presented in this section 
are the main result of this work. The essential fea- 
tures of the effective formalism developed here are 
that it is nonlinear and that it is not tied up in any 



way to a gradient expansion. Note that the nonequi- 
librium tensor 7 pi/ , from which the shear tensor IP" 
is obtained a posteriori as a quadratic function, sat- 
isfies a differential equation (|104[) instead of being 
an algebraic function of velocity gradients (as it is 
IP" in fluid dynamics) . In the colorless case, it was 
shown in [8l|, [82| that this feature results in a faster 
isotropization of the pressure as compared to second 
order hydrodynamics. Moreover, as opposed to the 
case of hydrodynamics, in the (colorless) effective 
theory the longitudinal pressure is positive through- 
out the entire evolution. We note that similar results 
were obtained in |77| within the so-called anisotropic 
hydrodynamics approach. We expect that these two 
results hold also in the present case including color 
degrees of freedom, although numerical simulations 
are needed to verify this. 

The color fugacities enter nontrivially in both 
sides of the total stress-energy tensor conservation 
equation, Eq. (|92l) . In the left-hand side they en- 
ter through the expression for the matter stress- 
energy tensor, Eq. (|9ip . while in the right-hand side 
they enter through the coupling to nonabelian fields. 
Moreover, both the hydrodynamic variables (u p , p) 
and the nonhydrodynamic variable 7 P " couple to the 
color fields and to £f m the evolution equation (11041) . 

From Eq. (|104|) . it is seen that the dynamics of 
the system is highly nontrivial. In part, this is be- 
cause the color chemical potentials must adjust to 
the flow as well as to the evolving gauge fields in 
order to make the system globally colorless (recall 
we obtained the evolution equation from the singlet 
sector of the kinetic equation). We emphasize that a 
nonvanishing color current does not imply that the 
system as a whole carries a finite color charge, be- 
cause the space-time dependence of the color chemi- 
cal potentials can be such that the total color charge 
vanishes [27| . 

The information about the constituents of the mi- 
croscopic theory is encoded in the transport coeffi- 
cients of the effective theory, and should in principle 
be computed from the former. However, the trans- 
port coefficients can also be treated as adjustable 
parameters. A well-known example is the case of 
fluid dynamics, which is usually derived from kinetic 
theory in the weakly coupled limit and under the re- 
laxation time approximation (to first order in r), 
but then used to describe strongly coupled matter 
@, [H • This is done by replacing the transport coef- 
ficients of the kinetic theory by those corresponding 
to strong coupling, which must be computed by dif- 
ferent means, for example, by using the AdS/CFT 
correspondence [H, |99|. A similar route can be 
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taken with the various transport coefficients that 
arise in our formalism (at linear order in j tJ these 
are (77,7-^,01,02)). In this regard, it is important to 
emphasize that the effective theory presented here 
is consistent by itself, independently of its deriva- 
tion from kinetic theory (this is also true in the case 
of second order hydrodynamics [H, [9{|), because it 
satisfies the Second Law and it is expected to be 
causal (the colorless version was shown to be causal 
in 0). 

The reason for considering the colorless part of the 
transport equation is that, as discussed in [32|, H(|, 
the color currents can persist in the plasma signif- 
icantly longer than the color charge density, which 
is neutralized rather fast. Therefore, a reasonable 
hypothesis is to assume that the dynamics of the 
system is determined by the singlet part of the dis- 
tribution function. What this means is that, even 
though the particles carry color (and thus interact 
with nonabelian fields), what we are actually de- 
scribing with the effective theory is the collective 
(or macroscopic) behavior of these particles, and it 
is this collective flow that is colorless. This hypothe- 
sis is physically well-motivated because one does not 
expect the plasma to be globally colorful, and due 
to this fact it has been adopted in previous studies 
dealing with chromohydrodynamics [26l-l29| and ki- 
netic theory [H, HH, [52| (see also |55|| for a related 
discussion of this issue in the context of the so-called 
"color hierarchy" transport equations). 

3. Yang-Mills equation 

In order to obtain a selfconsistent system of equa- 
tions for the variables (p, u^, , 7**", A afJ/ ) : the con- 
servation equations (1U21 and (|M|) and the evolution 
equation (|104[) must be supplemented with Yang- 
Mills equation for the gauge fields. For complete- 
ness, we write it down 

= nu»(cX + \d b ad CiCt) 

(107) 

IV. RELATION TO MATRIX 
CHROMOHYDRODYNAMICS 

Although our approach involving nonhydrody- 
namic variables is not tied up to a gradient expan- 
sion, for a local equilibrium state \^ vanishes and 



so the effective theory reduces to ideal fluid dynam- 
ics. The ideal fluid chromohydrodynamic approach 
has been discussed in [13, H3, Hll , and later on ap- 
plied to diverse studies mostly related to plasma in- 
stabilities and medium-jet interaction in the context 
of heavy ion collisions @, [H [H, H| . 

In this section we will compare the linearized ver- 
sions of our approach to linearized matrix chromo- 
hydrodynamics of Ref . [27j • We remark that beyond 
the linear order it is not possible to establish a simple 
mapping between these two approaches. In the ma- 
trix chromohydrodynamic approach, the standard 
procedure is to linearize the matrix equations 

D 2 ,T M " = 3 v ¥ vli (108) 

D M J M = (109) 

together with Yang-Mills equation in matrix fluc- 
tuations <5u M , Sp, 8J£, 5F avfl with respect to a given 
background. Aditionally, a relation between fluctu- 
ations of the matrix energy density p and the matrix 
pressure p is used. To the best of our knowledge, in 
previous studies based on the matrix approach the 
relation used in all cases is (Sp) a = c 2 s (8p) a . It is 
worth noting that the use of Sp — c 2 s 5p in our ap- 
proach is completely equivalent to the one adopted 
in the matrix formalism, as will be shown shortly. 
In studies dealing with chromohydrodynamics, usu- 
ally the ideal fluid case is considered; for a recent 
extension to the Navier-Stokes case see [26l . [29| . 



A. Generalities 

To understand the connection between ideal or 
Navier-Stokes chromohydrodynamics and our effec- 
tive theory, we first note that our development of the 
effective theory was ultimately based on Eq. (|10l) . 
which determines the coupling of matter to fields. 
In Eq. (fT0|) , the gauge fields are coupled to the trace 
of the matter energy-momentum tensor. Therefore, 
we end up with evolution equations for tr(T M "), and 
not for T M1/ itself, coupled to the Yang-Mills equa- 
tion. Instead, the chromohydrodynamics of [27J is 
based on T^, which requires the introduction of 
color matrices u M , p. 

The color current, defined in Eq. (|13[) . is then 
written as J M = nu M . In our formalism, the quan- 
tity n( a can be interpreted (at least at linear order) 
as the average color charge of a stream of colored 
classical particles. This interpretation for £ a can be 
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seen from its linearized equation of motion given in 
Eq. © 



.({+gnu»C J ab AlC 1 =0 



f 



(110) 



We see that it is identical to Wong's equation [57| 
(see also [IH [H, HH, [62j]) for an average classical 
color charge Q a = nCf, with the time derivative of 
the particle's trajectory replaced by u' 1 , i.e., the flow 
velocity. 

Having the interpretation for £ Q described above 
in mind, the linearized colored fluctuations of the ve- 
locity and the four-flow n of the matrix approach can 
be written in terms of the scalar fluctuations used 
in our approach as Su 1 ^ = (i a 5u^ and Sn a = nC,\ a , 
where £i a stands for the background value, which 
must be nonzero to avoid ending up in the case of a 
truly colorless system, as opposed to one which has 
color degrees of freedom but is in a colorless equilib- 
rium state. This mapping will be used in the next 
section to obtain the polarization tensor of the col- 
ored plasma. 



B. A simple example: The polarization tensor 

The polarization tensor characterizes the linear 
resp onse of the system to exte rnal perturbations 
[Ml45lM[^[^.[6l[6l [l^ . [l03l . and it is therefore 
an interesting quantity to compute in the formalism 
presented here. With the mapping between the ma- 
trix and our approach described in the previous sec- 
tion, it is straightforward to show that indeed the 
linearized equations of motion have the same struc- 
ture, and thus a very similar polarization tensor is 
obtained. However, the linearized equations are the 



same only if we set 7^" = or 7^" = —r\a^ v , which 
where the only cases studied so far |27H29j . 

Our formalism naturally incorporates higher order 
velocity gradients, since the dissipative tensor 
evolves according to a differential equation rather 
than being expressed as an algebraic function of ve- 



-7/(7 



jil' 



locity gradients, e.g. at first order as 7 
It is therefore interesting to investigate the role of 
higher order terms on the polarization tensor. Such 
terms are important, for instance, in the earliest 
stage of evolution of the QGP created in heavy ion 
collisions or in sit uations w here a fast parton goes 
through the QGP [l00l fl0ll|. 

In order to quantify the impact of higher order 
viscous terms on the linear response of the system 
to an external gauge field, we compute the polariza- 
tion tensor explicitly in our setting. The polarization 



tensor is defined through 

6J£ = - 



1 ah A bv 



(111) 



where is a small external perturbation. For 
simplicity, we shall consider an homogeneous, sta- 
tionary and colorless background described by n, 
and p. Using the mapping described above, 
namely 5u^ = (ladu 1 * and Sn a = n-Cia, together 
with (5p) a = CiaSp and 7^" = Cia7 Miy , our linearized 
equations become (covariant derivatives become or- 
dinary derivatives at this order) 

nd^Su^ + u^d^Sria = 

U^Spa + (1 + cDpd^uZ = 

c 2 s A^dJ Pa + (l + c 2 s )pu a d a 5uZ 

- nu a Fr + & al r = 

1 



(112) 



d v F» v = u»6n B 



-{wl P +l S a P ) 

+ n5u^ 



where we have put A 1 * 1 ' = g^ v + u^u v . For reasons 
that will become clear soon, in the above equations 
we use a generic squared speed of sound c 2 instead 
of the conformal value c 2 = 1/3, i.e. the relation 
between pressure and energy density perturbations 
reads Sp — c 2 Sp. 

Performing a Fourier transformation we can ex- 
press 5J% in terms of the gauge field perturbation 
and thus find the polarization tensor (the calcula- 
tion is very similar to that presented in Ref. 29]). 
We get 



1 ab 



>ab 



V 



(l + W 2 W A ){k-u) 2 
k 2 rt l u u - (fc • u) V 



+ (Wi + W 3 )[k 2 W 5 - k^k"(k ■ uf - k A u»u v 

(113) 

where 

Wi = -( k 2 + (c- 2 - l)(k ■ u) 



W 2 

w 3 



V 



{l + c? s )p(k ■ u)[l + iT„(k ■ u)] 

W 2 (l+4W 1 W i ) 
~3 + 3c^+4WW 4 



(114) 



W 4 = k 2 - (k ■ u) 2 
Wr = {k-u){u^k v + k^u v ) 
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and 



(l + cf)p 



(115) 



is the plasma frequency. Note that T^Y is diagonal 
in color space, as expected since as shown by Eqs. 
(|112| there is no mixing of colors at linear order, and 
transverse with respect to k^ . 

The result for that we obtain is the same as 
that obtained in Ref. (2(| but with the shear viscos- 
ity rj replaced by an effective one 



n 



1 + iT^(k ■ u) 



(116) 



The appearance of r) c s in place of r\ is quite natu- 
ral since 7V is precisely the relaxation time of the 
shear tensor IP" towards its Navier-Stokes value. 
We emphasize that this similarity with the chromo- 
hydrodynamic result holds only when linearizing the 
equations. The fully nonlinear evolution equations 
of the effective theory developed here, which do not 
involve hydrodynamic gradients, are different from 
those of Navier-Stokes chromohydrodynamics. 

To quantify the effect of higher order viscous terms 
on the linear response of the plasma to an external 
gauge field, we will work with the longitudinal part 
ti of the dielectric tensor e IJ . Similar analysis can 
be performed for the transverse part of the dielectric 
tensor, ex, but for brevity we shall only consider e^. 
We have (we suppress color indices and put fc M = 
(w,k) in what follows) 



value computed from kinetic theory of a Boltzmann 
gas (without color), which is given by 7v = 6rj/(sT). 
As a typical value for the temperature we shall use 
T = 200 MeV. 

The comparison made between ideal chromohy- 
drodynamics and kinetic theory carried out in Ref. 
(33j in the context of jet-induced instabilities shows 
that in order to achieve reasonable agreement be- 
tween both descriptions an effective speed of sound 
must be used in the former. To be consistent with 
previous studies we will show numerical results ob- 
tained with this effective speed of sound (instead of 



1/3), which is given by 



C * = 3 



2y \l + y 



\ (121) 



where y — V^ 2 /co. We note that this expression for 
c 2 is chosen to make the longitudinal dielectric func- 
tion obtained from the ideal chromohydrodynamic 
approach of Ref. [33[ identical to that obtained from 
kinetic theory in the leading-order HTL approxima- 
tion d [H Specifically, in the soft limit 
ui,Vk^ <C T and putting k = (fc, 0, 0) for simplic- 
ity we have [29| 



Hi 
k 2 



12a; pl 7?cff^ 

fc 2 sT 



LO 

2k 



In 



co ■ 



1 ~ i'° 



to — k 
lo + k 



in0(k 2 -lo 2 ) 



LO 



- k 



4fc 2 





LO + k 




In 








lo — k 





1 2 



1 



uo A 



and 



k 2 



In the rest frame — (1, 0, 0, 0) we get 

, ,2 



e L (uj, fc) = 1 \ 



1 



lo 2 V 1 - W 2 k 2 



(117) 



(118) 



where we can use that (1 + c 2 )p 
tropy density) to rewrite W2 as 

V 



w 2 



sTlo(1 + IT^LO) 



(1 - {Wx + W 3 )) 
(119) 

= sT (s is the en- 



(120) 



For illustrative purposes, and following |29j, we 
focus on the soft modes lo, Viz 2 <§; T and set -s/k 2 = 
2w p i and T = lOwpi. The relaxation time is set to its 



4fc 2 



ir 2 d(k 2 -io 2 ) 



+ in 



co 

2k 2 



In 



lo + k 



6{k 2 



(122) 



where 8 — > + and we used that 



In 



lo + k + iS 



lo — k + i5 



= ln 



- mQ{k 2 - lo 2 ) (123) 



with 8(x) the Heavyside step function. This shows 
that, in the soft limit, the modes with co > fc are 
undamped if 7v = 0. If the relaxation time does not 
vanish, the imaginary part of is nonzero even for 
lo > k. We will now show some illustrative examples 
of this feature. 

Figure [T] shows the real and imaginary parts of cl 
as a function of co /k for r\j s = 0.3. This value for 77/s 
is on the high side in terms of fitting viscous fluid 
dynamics results to RHIC and LHC data It 
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is seen that there arc significant differences between 
the longitudinal dielectric function computed with 
both values of the relaxation time. The most no- 
ticeable effects are seen on the imaginary part of cl , 
which, for lj < k, is smaller in the case with nonva- 
nishing r v . This indicates that, in this range of fre- 
quencies, the induced color excitations decay more 
slowly as compared to the case with t w — 0. This 
behavior can be understood by recalling the physi- 
cal meaning of t v as the relaxation time of the shear 
tensor IP 1 ' towards its Navier-Stokes value —r}G^ v . 
If the value of r w is increased, then hydrodynamic 
fluctuations will decay more slowly. Since hydrody- 
namic fluctuations are coupled to color fluctuations, 
the latter will decay more slowly as well. This result 
is in agreement with those of Ref. [5l| obtained from 
kinetic theory with a BGK collision kernel, showing 
that the addition of hard-particle collisions slows the 
rate of growth of QCD plasma unstable modes. 

As shown in Figured! the results that we obtain 
for u> > k show that the damping of color excita- 
tions in this frequency range is completely different 
according to whether is zero or not. As expected, 
if r n =0 there is no damping in this frequency range. 
This feature stems from the analytic structure of the 
longitudinal dielectric function in the regime where 
u,k < T, as given by Eq. (I122p . In this limit, 
the imaginary part of is proportional to the step 
function 9(k 2 — lj 2 ) (see e.g. |62[). On the con- 
trary, if tv 7^ then those color excitations with 
lj > k become considerably damped, with a damp- 
ing rate which falls off steeply with increasing fre- 
quency. Similar results were obtained in Ref. [49( 
for QED dispersion relations obtained from kinetic 
theory with a BGK collision term. There it was 
found that when collisions are included, the longitu- 
dinal dispersion intersects the light cone lj = fc, in 
contrast to the case of collisionless dispersion where 
lj > k for all k. In a collisionless plasma, Lan- 
dau damping, which is possible only for lj < k, 
is the only damping mechanism, and thus plasma 
waves are undamped. In contrast, collisions intro- 
duce an additional damping mechanism for plasma 
waves (see e.g. [5l|). We emphasize that a first 
order hydrodynamic formalism, in which the shear 
tensor relaxes instantaneously to its Navier-Stokes 
value, can not completely account for such damping 
of plasma waves. 

We now briefly discuss the influence of the value 
oirj/s on the longitudinal dielectric function. Figure 
[5] shows the real and imaginary parts of for a 
smaller value of the viscosity-to-entropy ratio than 
before, namely r//s = 0.15. It is seen that the impact 



of a nonvanishing relaxation time on et decreases 
with decreasing values of r\j s. For the value r\j s = 
0.15, the effect of tv on is still significant, and we 
still see that excitations with u > k are damped due 
to collisions. Although not shown, we find that for 
r)/s < 0.08 the difference between the longitudinal 
dielectric function obtained with a vanishing or a 
nonvanishing value of is hardly appreciable. 
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FIG. 1: (Color online) Real (upper panel) and imaginary 
(lower panel) parts of the longitudinal dielectric function 
e_L as a function of frequency, for rj/ s = 0.3. The values 
of the parameters are set to k = 2ca p i, T = 10o;pi, tv = 
6v/(sT) or t w = 0, and T = 200 MeV. The results are 
obtained using the effective speed of sound given by Eq. 

ifmj) . 

The effective formalism discussed in this work is 
phenomenological and involves several approxima- 
tions. In spite of this, it is interesting to qualita- 
tively discuss possible implications of our results for 
the phenomenon of jet quenching in heavy ion col- 
lisions. For this, and considering the already dis- 
cussed limitations, we take the view that our formal- 
ism constitutes an appropriate model to understand 
some features of the response of the QGP to a fast 
moving parton that crosses it. 

There are two main energy loss mechanisms which 
contribute to energy loss: radiation of soft glu- 
ons and collisions involving the exchange of hard 
(~ T or larger) and soft (~ 2u p i) momenta (13 , 
Il4j . The dominant source of energy loss is gluon 
bremsstrahlung, although the contribution of colli- 
sions to the total energy loss is significant, particu- 
larly when attempting to fit the results of theoretical 
models of jet quenching to data [HI, [2 EH, H^]. In 
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FIG. 2: (Color online) Real (upper panel) and imaginary 
(lower panel) parts of the longitudinal dielectric function 
e_L as a function of frequency, for rj/s = 0.15. The values 
of the parameters are set to k = 2uu p \, T = 10aj p i, tv = 
6j7/(sT) or -rv = 0, and T = 200 MeV. The results are 
obtained using the effective speed of sound given by Eq. 
llSIIl- 



connection to our results, we note that the contri- 
bution of soft collisions to the energy loss can be 
dire ctly calcu lated from and £t (see for example 
mElSEol]). We sha11 not discuss in detail this is- 
sue here, but just mention that a smaller imaginary 
part of £l will result in a decrease in the energy loss. 
Our results then show that a nonvanishing value of 
the relaxation time t„ is expected to lead to a size- 
able reduction in the energy loss. 

As a final remark, we note that the radiation spec- 
trum and hence the energy loss of a hard parton 
crossing the QGP is also modified by the dielectric 
polarization of the medium (this is known as the Ter- 
Mikaelian effect jl05l | - s ee |l06l| for an extension to 
QCD). As noted recently jl04| . the effect of radiation 
damping occurring in an absorptive medium on the 
spectrum of radiated gluons is particularly interest- 
ing and might lead to sizeable effects on energy loss 
related phenomena. The polarization tensor derived 
from the formalism presented here naturally incor- 
porates damping. A detailed study of the influence 
of tv on the energy loss of fast partons for the condi- 
tions prevailing in heavy ion collisions at RHIC and 
LHC is left for future work. We emphasize, however, 
that the richness of the effective formalism presented 
here lies in its nonlinear character, which is not re- 
flected in the polarization tensor but may become 



relevant when dealing with parton energy loss phe- 
nomena . 



V. SUMMARY AND OUTLOOK 

In this work we have obtained from the kinetic 
theory of nonabelian plasmas an effective model de- 
scribing the evolution of a system composed of col- 
ored particles interacting with nonabelian classical 
gauge fields. The link between the one-particle dis- 
tribution function of colored particles in the kinetic 
description and the variables of the effective the- 
ory is determined by the entropy production varia- 
tional method. The closure provided by this method 
does not rely in any way on a gradient expansion in 
macroscopic variables and can therefore be applied 
even when these gradients are large. 

In order to compare the developed effective the- 
ory with chromohydrodynamic formalisms based on 
the usual gradient expansion, we have calculated the 
longitudinal dielectric function of the plasma. Us- 
ing typical values of the plasma parameters appro- 
priate for the QGP, together with an effective speed 
of sound chosen to reproduce the longitudinal dielec- 
tric function of hard-thermal loop kinetic theory, we 
have found that the relaxation time for the shear ten- 
sor has a strong influence on the dynamics of color 
fluctuations, in agreement with the results of kinetic 
theory including collisions among the hard partons. 
The implications of such changes on the evolution 
of color excitations on phenomena relevant to heavy 
ion collisions, particularly on jet quenching, deserve 
further investigation. 

The formalism presented here is a simplified model 
of the dynamics of color fields during the early and 
intermediate stages of heavy ion collisions. It may 
be useful to shed light on issues which would require 
intensive simulations in a microscopic approach, for 
example the magnitude of the back reaction of the 
particles' flow on the gauge fields. If the color fields 
eventually die out, the effective theory goes over to 
second order fluid dynamics if the velocity gradi- 
ents are small, so that the effective theory could be 
used to describe (starting from suitable initial condi- 
tions) the evolution of the fireball created in a heavy 
ion collision from very early times (> 0.2 fm/c) till 
freeze-out in a unified, albeit simplified, way. 

Concerning the dynamics of color fields at early 
and intermediate times, our formalism could be used 
to study plasma instabilities and its effect on the 
evolution of matter created in heavy ion collisions. 
It would be particularly interesting to solve numer- 
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ically the full nonlinear equations of the effective 
theory presented here for the conditions prevailing 
in heavy ion collisions, and to compare the results 
to those obtained by a microscopic approach. To 
carry out this program, the inclusion of hard gluons 
into the model is certainly required for a realistic de- 
scription of the physical processes involved at those 
stages. Work is in progress along these lines. 

At this point we would like to comment on the 
limitations of our approach in connection to possi- 
ble applications to describe the early-time dynamics 
of color fields in heavy ion collisions (some of the 
issues discussed here are also relevant for parton en- 
ergy loss). We think that it is clearer to distinguish 
between limitations inherent to our approach (that 
are either truly unsurmountable or else very difficult 
to address) and simplifying hypothesis that could be 
relaxed in the future. 

We start by discussing those limitations that are 
inherent to our approach. 

The effective theory presented here constitutes a 
simplified model of the true dynamics given by ki- 
netic theory with a linear collision term, and there- 
fore the equations derived in this paper should be 
applicable for similar time scales. The kinetic the- 
ory description of the early stage of heavy ion colli- 
sions is valid for times > Qj 1 (~ 0.2 fm/c at RHIC) 
when particles having transverse momenta greater 
than Q s are formed out of the color fields |{§ llOj . 
However, there are some limiting factors that should 
be considered. 

The first one was already hinted to and refers to 
the use of a linear collision term. Although one 
could, in principle, write down the variational equa- 
tions of the EPP for a transport equation with non- 
linear collision kernels, there is no much prospect 
of being able to solve them. However, this may 
not be a serious issue, because several studies have 
shown that the Boltzmann-Vlasov equation coupled 
to Yang-Mills equation provides a fairly reliable de- 
scription of the early stag es of a heavy ion colli- 
sion 0, d m Mm, MM EM M HI. The 

short-range interaction between hard partons has 
been incorporated in the kinetic models only re- 
cently [H, Hi M m H3 . 

The second point is that the range of applicability 
of the effective theory developed here is not deter- 
mined by the magnitude of velocity gradients, but 
rather by whether the dynamics of the system as 
given by (linearized) kinetic theory can or can not 
be described by few variables (including nonhydro- 
dynamic ones) coupled to classical gauge fields. As it 
happens with other approaches to the closure prob- 
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is difficult to precisely establish a priori the range 
of validity of the resulting effective theory. We be- 
lieve that, as it happens with fluid dynamics, our 
formalism may prove useful to describe some stages 
of a heavy ion collision provided the transport coef- 
ficients are suitably chosen. In any case, this point 
should be settled by comparing the results obtained 
from the effective theory to those obtained from ki- 
netic theory 

As an issue that can be improved in the future, we 
mention first the fact that here we deal with excita- 
tions of scalar fields (and not spinors), and second 
that we do not take into account the hard gluons. 
In the early stage of heavy ion collisions, the hard 
gluons are, as the hard quarks, coupled to the soft 
gluons and therefore should also be described by a ki- 
netic equation with its corresponding collision term. 
This term would involve not only interactions among 
hard gluons themselves but also between hard glu- 
ons and the excitations of the scalar fields. The lat- 
ter coupling should also be reflected in I co / of Eq. 
(fTTj) . The inclusion of hard gluons interacting with 
the classical gauge fields is mandatory for the model 
to be applicable to the early stage of a heavy ion 
collision, in which strong fields decay into gluons to 
eventually form the QGP. 

Another point that can be addressed in the future 
is the possibility of performing a quadratic expansion 
around an anisotropic (in momentum space) distri- 
bution function instead of the isotropic distribution 
function we use here. This idea was exploited in 
[77l . [78j to obtain an effective theory that can handle 
the very large anisotropies in momentum space that 
are present at early times in heavy ion collisions, 
and moreover reproduces both the ideal fluid and 
the free-streaming limits. The correct description 
of both regimes (which are ultimately determined 
by the value of the Knundsen number) may be a 
relevant issue when dealing with plasma instabili- 
ties. Although our formalism can handle very large 
anisotropies as well, at present it is not clear how 
well can it describe plasma instabilities, so further 
studies are certainly needed to settle this point. 

Finally, we plan to include stochastic terms into 
the evolution equations of the effective theory jl07l . 
Il08j . This would allow us to address in a sim- 
ple model setup the important question on the fate 
of fluctuations and their impact on observables at 
RHIC and LHC. 
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